o 
o 
o 






NYU-99-06-02 

Violation of the Greisen-Zatsepin-Kuzmin cutoff: A tempest in a 

(magnetic) teapot? 
Why cosmic ray energies above 10 20 eV may not require new 

physics 



Glennys R. Farrar and Tsvi Piran* 

Department of Physics, New York University, NY, NY 10003, USA 
* Permanent address: Racah Institute of Physics, Hebrew University, Jerusalem, Israel 



<N ■ (February 1, 2008) 

m 

vo, 

Q ! Abstract 

Pu,; 

O ' The apparent lack of suitable astro-physical sources for the observed high- 

est energy cosmic rays (UHECRs) within ~ 20 Mpc is the "GZK Paradox". 
We constrain representative models of the extra-galactic magnetic field struc- 

X' 

$— i ' ture by Faraday Rotation measurements; limits are at the //G level rather 

than the nG level usually assumed. In such fields, even the highest energy 
cosmic rays experience large deflections. This allows nearby AGNs (possi- 
bly quiet today) or GRBs to be the source of ultra-high energy cosmic rays 

without contradicting the GZK distance limit. 
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In 1966, Greisen, Zatsepin and Kuzmin JT| (GZK) pointed out that high enough energy 
protons degrade in energy over cosmologically small distances due to photopion production 
from the cosmic microwave background. Less than 20% of protons survive with an energy 
above 3 x 10 20 (1 x 10 20 ) eV for a distance of 18 (60) Mpc 0. Ultra-high energy (UHE) 
nuclei and photons lose energy even more readily. Yet more than 10 cosmic rays (CRs) 
have been observed with nominal energies at or above 10 20 ± 30% eV |||4]] with the Fly's 
Eye event having 3.2 x 10 20 eV [g|- Elbert and Sommers, using search criteria based on the 
Hillas condition and other reasonable expectations as to the properties of sources capable 
of accelerating protons to these high energies, found no sources within 50 Mpc of Earth ||. 
Thus, it is widely believed that the GZK limit must somehow be violated to explain the 
origin of the observed UHECRs. 

The need for GZK violation is based on the assumption that UHECRs and photons from 
the same source should arrive from the same general direction and with only a moderate 
arrival time difference. If this assumption is not valid, the observed UHECRs could be pro- 
duced by sources within the GZK distance that are not in the direction of the UHECRs and 
moreover may have evolved significantly between the UHECR production and the emission 
of the photons we now observe. In this case, the existance of particles with ultra-high ener- 
gies would not raise a paradox if one or more candidates exist for a "waned" source, within 
the GZK distance. 

The coincidence of an observed UHECR and its astrophysical source depend on the 
angular deflection of the UHECR being small. The UHECR is deflected by the magnetic 
fields it encounters which are generally assumed to be of order nG, leading to expected 
deflections of order a few degrees or less: 



59 ~ 0.5 o yjD Mpc \ Upc (B nG /E 20 ) (1) 

being the rms deflection of a proton of energy .E20IO 20 eV traveling a distance D through 
randomly-oriented patches of magnetic field having rms value B and a scale length A ||. 
The corresponding difference between the arrival times of photons and UHECR's, r arr ~ 



(56) 2 D/2c, is small on astrophysical scales, ~ 10 4 yr taking D = Lqzk ~ 20 Mpc for 
a E 2 q = 3 proton in a nG average field. With these values both angular and temporal 
correlations are expected between the photons and UHECRs of a given source. 

However, as we show here, extra-galactic magnetic fields (EGMF)'s are plausibly larger 
by several orders of magnitude than is generally assumed, in which case the angular deflection 
of UHECRs and the time delay between a UHECR and photons from the same source have 
been grossly underestimated. With larger fields there is no angular and temporal correlation 
between a UHECR and its source. The long time delay between the arrival of UHECRs 
and of photons implies that the source of the UHECRS could be unremarkable or even 
undetectable today. The expected cutoff in the UHECR spectrum is modified as discussed 
below. 

Several groups have recently explored the possibility that the EGMF has a large scale 
structure akin to that observed in simulations of dark matter, assuming mean fields in the 
0.1-1 /iG range |7|||. Results of our analysis are consistent with these detailed numerical 
simulations, and clarify the generality of a large EGMF, >few tenths yuG, relieving the 
concern that those results are artifacts of the particular spectra of magnetic inhomogeneities 
employed (e.g., Kolmogorov [0, log- normal []8|]). Subsequent work P,|T0"[ explored these field 
ansatze to find what parameters give the best fit to the observed spectrum and angular 
distribution of CRs above 10 19 eV. We instead focus on the puzzling highest energy events, 
those above 10 20 eV. We show that the broad angular distribution and absence of identifiable 
sources within the GZK volume is not a problem for either GRB or AGN sources. Our 
discussion focusses on the interpretation of UHECRs as protons since for nuclei the energy 
attenuation distance via photodissociation is shorter than the GZK length for protons of the 
same energy. However our results can be applied to nuclei as well. 

The structure of this letter is the following. We first determine the Faraday rotation 
limits on the EGMFs. Next we examine the trajectories of UHE protons in these EGMFs 
and debunk the notion that the deflection angles of 10 20 eV protons must be of order a few 
degrees or less. This removes the necessity of an angular and temporal correlation between 



UHECRs and their astrophysical sources. Finally we show that AGNs or GRBs active in 
the local supercluster in the past 10-100 million years can account for the observed UHECR 
flux. We conclude with a summary and some observational tests of these scenarios. 

EXTRAGALACTIC MAGNETIC FIELD: The prime constraint on the EGMF arises 
from the Faraday rotation of light from distant quasars JTTJ , for which cosmological effects 



must be included. The rotation measure of a source at a redshift z is 

where ft.75 is the Hubble constant in units of 75km/sec/Mpc, distances are in Mpc, and 7105 
is the electron density in units of (fib /0. 05) (3Hq/87iG) /m p ; fi&, m ,A are the ratios of bary- 
onic, matter, and vacuum energy densities to the closure density. The comoving distance 
increment, dr c , is related to the physical distance by dl = dr c /(l + z'). Additional redshift 
factors arise due to the RM's quadratic dependence on frequency and the possible cosmo- 
logical evolution of the electron number density n e (z) = 71-05(1 + z) Pe and the magnetic field 
strength B(z) = Bq(1 + z) PB . The cosmological "dimensionless effective distance" d p (z), 
normalized so that do(z) measures the co- moving distance in terms of the Horizon distance 
2c/ ' Hq and with p = p e + p B — 3, contains the net effect of the cosmological evolution: 

*w - f & ' (1 + zr/2 r . ( 3 ) 

In the models below, for z = 2.5, Q\ = 1 — Q m and Q m = 1 (0.3), we encounter d 2 — 
1.85(2.98), d = 0.46(0.68), and cL 3 = 0.14(0.18). 

Faraday rotation measurements for constraining the EGMF make use of quasars 
with (z) m 2.5 and yield RM<5radm~ 2 [11]. Eqn (0) can be written RM5 = 
400 n 05 B\\^ G h 75 d P (z), where RM(z) = 5 RM 5 rad m -2 . We now consider two extreme 
models of the structure of the EGMF. 

Randomly Oriented Patches: A commonplace model, which leads to eqn ([I]), posits that 
the EGMF consists of domains of constant but randomly oriented field, much like in a 
ferromagnet, with present rms strength Bq and characteristic size A | ll| . Assuming that n e 



and A are constant in co-moving coordinates, and that the energy density of the magnetic 
field scales like radiation, implies 



R HA n RM5 l V^ (A\ 

W5 n 05 \l A Mpc a ^\ z > 
where we have taken the mean number of patches in the path to the quasar to be r c (z)/X. 
Assuming n e « rij, and baryonic closure of the Universe, i.e., n 05 « 20, implies nG fields 
for A ~ 1 Mpc |iT] and hence the small deflections generally assumed in discussions of 
UHECR trajectories. However n e should be taken at least a factor of 20 smaller than in 
those estimates, since f2& < 0.05. Moreover Qj, has contributions from neutrons, and only 
electrons in ionized gas are relevant to Faraday rotation, so this further reduces n e . Taking 
the more realistic value n 05 « 0.3 gives B in this model of order 0.5 (0.4) RM 5 (0.3/ra 5)yuG, 
for the flat Q m = 1 (0.3) cosmologies. 

Sheets and Voids: It is more likely that the scale structure of the EGMF consists of 
randomly oriented sheets of field, presumably associated with the sheet-like and filamentary 
concentrations observed in the matter distribution, separated by relative voids. Idealize this 
to sheets consisting of layers of thickness L$ within which the field has a constant magnitude 
Bq and random orientation, separated by voids of thickness Lo(~ 50 Mpc), with Ls <C Lq. 
Observations of high redshift clusters provide strong evidence that Q m is small and the sheets 



and voids of matter were largely in place at z ~ 2.5 | I2| . Therefore we expect no scaling 
of B, n e , and Ls with redshift in the sheet-void model, for the redshifts < 2.5 relevant for 
Faraday rotation observations. Thus eqn (ffl) can be used to obtain B in this model by 



replacing (^(z) — > d-s(z) and multiplying by JLq/Ls- Note that Bo and n e stand for the 
present rms field strength and electron density within the sheets and not averaged over the 
voids, so n e can be expected to be enhanced compared to the random-patches model by the 
factor Lq/ Ls- Therefore in the sheet-void model we find B of order <i2(2.5)/(o?_3(2.5)v / 50) 
times that of the random patches model. This is a factor w 0.9 in both cosmologies. If the 
sheets themselves consist of many patches of randomly oriented field of typical size A, B 



should be reduced by a factor JX/L 



The Faraday rotation estimate, B »s 0.5 RM 5 /xG, is only an upper limit. However, this 
value is comparable to the ~ 1.5/xG field argued to exist at the core of the local (Virgo) 
Supercluster [13j and to field intensities of ~ 0.2 /xG observed in Abell 2319 and Coma (see 



fTT| for a review and |14|], which appeared after our paper was submitted, for further evidence 
for such high magnetic fields.). A /xG field corresponds approximately to equipartition 
between the magnetic energy and the gravitational and thermal energies of the supercluster. 
The consistency of the inferred field from these disparate approaches increases our confidence 
in the conclusion: The EGMF, at least in the local supercluster, is plausibly of order a 
few tenths /xG rather than nG as has been traditionally assumed. In the random patch 
approximation the energy density of this field is comparable to that of the CMBR. Therefore 
we focus primarily on the more likely "sheets and voids" model. 

TRAJECTORIES OF UHE PARTICLES: The Larmor radius of a proton with an energy 
i?2ol0 20 eV in a constant orthogonal field, B±, is 

R L = 0AlMpcE 20 /B ±llG , (5) 

corresponding to a deflection angle of order 59 « 0.5° B^ nG Am pc / 'E 2 o when travers- 
ing a distance A <C R^. In the random patches model, we can express the deflection 
angle directly in terms of the Faraday rotation measure allowing the patch-size and its 
uncertainty to be eliminated. For a source at a distance D ^> A the relationship is 



80 rms ~ 60° JDm P c R-M 5 /(ri 5 E 2 o h 7 ' 5 ), as long as this deflection is small. An analo- 
gous relationship can be obtained for the sheet- void model, however it is only applicable for 
cosmic rays which traverse many sheets, which is not the case for UHECRs since the GZK 
length is less than or of order the void size. 



Vallee ]H| quotes an average coherent enhancement >1.2 /xG in the central 10 Mpc 
region of Virgo, our local supercluster. There are surely some inhomogeneities at the 0.1 
- 1 Mpc scale, if only from the galaxies themselves. In such a field, the minimum energy 
above which UHECRs have approximately rectilinear motion over a distance of 10 Mpc 
can be estimated by requiring ([]]) to evaluate to < 10° with scale size A = 0.1 Mpc and 
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< B >= l.2/j,G. Such small deflection requires E>6 10 21 eV, far larger than in the field 
configurations considered in refs. P,[T0||, for which particles above (1 — 2) 10 20 eV display 



rectilinear motion. Given that fields may be larger than previously assumed, we should not 
rule out the possibility of magnetic confinement, diffusive motion, or large angle scattering 
for even the highest energy CRs thus far observed. 

IMPLICATIONS FOR UHECR SOURCES: Now we apply these results to possible 
sources of UHECRs. We consider either continuous sources (specifically AGNs) or bursting 
sources (specifically GRBs). In addition to being capable of accelerating the UHECR to the 
required energy the sources must satisfy two additional constraints. First, the effective num- 
ber of sources, N e ff, within the effective volume contributing to UHECRs on Earth, K//, 
must be > 1. If N e ff < 1, it is improbable that a source exists within the GZK distance 
during a time interval such that its cosmic rays arrive at Earth in the past 30 years. Second, 
the sources should be energetically capable of producing the observed flux of UHECRS on 
Earth. 

The effective number of active sources within the effective volume satisfies N e ff = V e ff x 
p/e s for constant sources, and N e ff = V e ff x Tr arr for bursting sources. Here p and r are 
the density and lifetime of active sources, r arr is the spread in UHECR arrival times, T is 
the rate per unit volume for bursting sources and e s is the "duty-factor" of a long-lived 
but not eternal source, e.g., e s = r arr /min(r, r arr ) for an AGN. The effective volume, K// 5 
depends on the magnetic field model. In the sheet-void model, the contributing region is 
the local supercluster within a distance of Lqzk'- if the cosmic ray experiences a few large 
angle scatterings as it travels from source to Earth, V e ff — > tiLq ZK D L sc, where D LS c ~ 10 
Mpc is the thickness of the local supercluster. 

The spectrum of UHECRs above 10 18 ' 5 eV can adequately be described by E 2 j(E) m 
3 eVcm _2 s _1 str" 1 , corresponding to an energy flux per decade of energy, $ b s , which is 
about 3 10 45 ergMpc _J yr _1 str _1 . The energy flux per decade in UHECRs, produced by N e ff 
sources is: 



$ = e CR E" e -^H, (6) 

4vr V eff 

where Eiq is the energy produced per energy decade in other forms (e.g. gamma-rays) by 
the source and ecu is the relative efficiency of producing an equal energy per decade in 
UHECRs. For a continuous source, E 10 = PioTqzk, where P w denotes the power per decade 
and t gzk = L GZK /c. 

AGNs and Hot Spots in Giant Radio Galaxies are the canonical example of sources that 
satisfy the Hillas acceleration conditions. However, there are no suitable AGNs within the 
GZK distance in the directions of the observed UHECRS. If large deflection angles are the 
norm, the alignment requirement is relaxed and all powerful radio sources within the GZK 
distance become candidates. Elbert and Sommers || list eight such high flux radio galaxies. 
Among them M87 (and possibly Cen A) are possibly strong enough to satisfy the Hillas 



condition [ I7| . These objects and others in the list might also have been more powerful 
AGNs in the past and thus are candidate sources @ in this picture. 

The UHECR flux from AGNs is computed using eqn || The power of a moderate 
AGN, in the energy decade corresponding to gamma rays, is of order 10 45 erg s _1 . N e ff in 
V e ff over the past GZK time must be > 1, as discussed above. The energy flux of UHE- 
CRs observed on Earth today due to a single AGN in our local supercluster is therefore 
$ ?s 10 48 erg Mpc~ 2 yr _1 €cr€ s . Hence from an energetics standpoint there is a comfort- 
able margin for inefficiency or inactivity in the UHECR production by AGNs in the local 
supercluster. 



If M87 is indeed powerful enough to accelerate UHECRs fll7 |, then it is possible that it 



produces all the observed UHECRs. The different arrival directions simply reflect different 
trajectories in the EGMF. An alternative scenario is that no source within the GZK limit (or 
only M87) is active today, but such sources have been active in the past. This implies that 
the UHECR-producing lifetime, tagn, must be <r arr ~ t gzk . If t agn ~^> t gzk , it would 
be unlikely that a source which was active recently enough to have produced an observed 
UHECR, would no longer appear active. On the other hand, if tagn "C t GZ k, the only way 
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to have N e ff > 1 is if pAGN^efj 3> 1- If this were the case, we should see many objects 
which might have been good sources during the last GZK time, while in fact we see just a 
few. For example, if tagn/tgzk = 0.1, the probability is ~ 43% that eight sources which 
have been active within tqzk would all be inactive today, and the probability is ~ 38% 
that 1 out of the eight would be active. Thus the likelihood of finding the current situation 
would be rather large. Therefore, we conclude that tagn ^ Tqzk ~ 6 10 7 yr. Remarkably, 
independent estimates of AGN lifetimes are in the 10 7 — 10 8 yr range [18], so the AGN source 
model survives a highly non-trivial requirement. 

GRBs: GRBs might be beamed with an opening angle 9grb ~ 1/10. With such an 
opening angle we observe only a fraction GRB /A of all bursts. Let T GRB be the overall rate of 
GRBs per unit volume. The observed rate of GRBs, {9 GRB /4)T GRB = T_g 1CT 9 yr -1 Mpc -3 , 
is independent of beaming and thus fairly well determined; if gamma ray bursts follow star 
formation T_g rs 1 |16||. Uniformity requires that the spread in arrival times of UHECRs 



from a single source is long compared to the time between contributing GRB events. If 
magnetic fields are of order nG or less, UHECRs experience small angular deflections. Only 
GRBs beaming toward Earth contribute, r arr rs 10 4 yr, and N e ff is marginal [IS . 



However if the magnetic field is larger and UHECRs experience large deflections, the 
beaming fraction suppression factor is removed. In addition, r arr increases because of the 
greater path length; it is replaced by the GZK attenuation time, tqzk- If the cosmic ray 
experiences a few large angle scatterings as it travels from source to Earth, N e ff fh 10 5 r_ 9 
for E = 3 10 20 eV, and the number of sources is adequate to assure our epoch is not unusual. 
Eqn | gives the energy flux per decade in UHECRs. Taking E 10 « 10 52 ergs, L GZ k ~ 60 Mpc 
(since E20 ~ 1 dominates the spectrum of UHECRs) and making the optimistic assumptions 
that [d GRB /A)~ 1 = 200 and e CR = 1, gives $ « 10 46 erg Mpc~ 2 yr _1 str _1 . Thus, Gamma Ray 
Bursts are viable, but not by a large factor, from this standpoint. 

PREDICTIONS OF THIS SCENARIO: Several clusters of 2-3 events having energies 
from rs 4 — 30 10 19 eV have been observed |J. If the ultra-high energy part of the cosmic ray 
spectrum consists of protons which experience large deflections due to few-tenth-/zG fields 
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in the local supercluster, this must be a statistical fluke. Even if UHECRs are produced 
in a handful of locations within the GZK distance, this initial localization will not survive 
subsequent deflection, which is large and strongly energy dependent without fine-tuning of 
the fields. Thus clustering or the possible directional correlation with distant compact radio 
quasars |HJ would be evidence against the picture advanced here. 



If M87 is the source of all UHECRs, several consequences follow from its being located 
near the center of the local supercluster, at a distance of 20 Mpc which is approximately 
the GZK distance for a 3 x 10 20 eV proton. When better statistics allow the spectrum to 
be investigated at higher energy a cutoff in energy should be evident in the data. Addi- 
tionally, if, as expected, the effective scattering distance is not too small compared with the 
supercluster dimensions we should find an "in-out" asymmetry with respect to the center 
of the local supercluster. Southern Hemisphere detectors, looking outward with respect to 
the local supercluster, would be expected to see a lower UHE flux than observed in the 
Northern Hemisphere. Since not much increase in pathlength can be tolerated for particles 
originating near the GZK distance, we could infer in this case that magnetic fields are not 
as strong as are allowed by Faraday Rotation and equipartition limits. If, instead, GRBs or 
multiple AGNs are the source of UHECRs, the energy cutoff would be softer than with a 
single source but the in-out asymmetry might still appear if the sources were concentrated 
near the supergalactic center. 

A striking aspect of this picture is that observed CRs come mainly from the local su- 
percluster: at ultra-high energy their range is limited by the GZK effect, while at lower 
energy they are confined near their sources by the relatively strong fields. This may lead to 
a unified explanation for the knee and ankle structures in the CR spectrum; a more complete 
treatment of the overall spectrum remains to be addressed in future work. 

In summary, we have shown that the highest energy cosmic rays can originate in the local 
supercluster, either from an AGN which now may appear past its prime or from gamma ray 
bursts which have occurred in the last 10-100 million years. We presented two distinct 
models of the structure of extragalactic magnetic fields, and recalled observational data on 
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the field in our own and nearby superclusters, which indicate that extragalatic magnetic 
fields can be of order /zG rather than the nG heretofore generally assumed. This radically 
changes the picture of ultra-high energy cosmic rays and implies that even the highest 
energy cosmic rays observed to date have strongly bent trajectories such that i) there is no 
directional correlation with the source; ii) there is no temporal correlation between UHECRs 
and photons from a given source, on time scales relevant for identifying the sources, even for 
AGNs, and iii) the power in UHECRs which can be supplied by the relevant superposition 
of sources contributing over a GZK time, ps 10 8 yr, is easily adequate for AGNs and may 
barely be adequate for GRBs. 

The research of GRF was supported in part by NSF-PHY-99-96173. We thank P. Bier- 
mann, C. Heiles, P. Kronberg, and T. Kolatt for valuable contributions. 
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